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Abstract 

We will prove the presence of chaotic motion in the Lorenz five-component 
atmospheric system model using the Melnikov function method developed 
by Holmes and Marsden for Hamiltonian systems on Lie Groups. 
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1 Introduction. 

The first model equations for the atmosphere are the so called primitive equa- 
tions (PE). This model allows wave- like motions on different time scales. The 
slow motions which have a period of order of days are called Rossby waves and 
the fast motions which have a period of hours are called gravity waves. The 
question of how to balance these two time scales lead Lorenz || to introduce a 
simplified version of the (PE) model, the so called five-component model. This 
is a system of five differential equations which couples the Rossby waves and 
gravity waves. This system turns out to have a Poisson formulation on R 5 first 
discovered by Bokhove EL We shall find a Poisson diffcomorphism between 
the Poisson structure of Bokhove and the product structure on se*(2) x M 2 , 
where the first factor is endowed with the Lie-Poisson structure and the second 
with the standard symplectic structure. Using this diffeomorphism, the five- 
component Lorenz model takes on a form appropriate for the application of the 
Melnikov function method on Lie groups developed by Holmes and Marsden Q . 
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2 The geometry of Lorenz simplified model of 
Rossby-gravity wave interaction. 

The model introduced by Lorenz in || is described by the following set of 
differential equations: 



Xi = 
X2 = 
X3 = 
X4 = 
x 5 = 



-X2X3 + £X 2 X 5 
XIX3 — 6X1X5 

~X\X 2 

-x 5 

Xi + SX1X2, 



(2.1) 



where the variables X4, x§ represent the fast gravity wave oscillations and Xx, 
X2, X3 are the slow Rossby wave oscillations, with the parameter e (related to 
the Rossby number) coupling the two sets of variables. 

In B, Bokhove wrote the system (|2.ip in the following Hamiltonian form: 



x = {x,H}x, 
where the Hamiltonian function is given by 



H(xx,x 2 ,x s ,X4,x 5 ) = ]^{x\ 
and the Poisson bracket {•, - }i is defined by 

{f,9h = 



2x 2 



X\ 



df_dg_ _ df dg 
8x2 dx3 8x3 8x2 
df_dg_ _ df dg 
8x3 dx\ dx\ dx3 
df_dg_ _ df dg 
dx§ 8x4 dxi dx^ 



-•''2 



df_dg_ _ df dg 

dx$ dx2 dx2 dx$ 
df_dg_ _ df dg 
dx\ dx§ dx$ dxi 



On the space se*(2) x R 2 consider the Poisson bracket given by the product 
bracket of the Lie-Poisson bracket on se*(2) and the Poisson bracket on R 2 
induced by the standard symplectic form on R 2 . Denoting by (hi, 112, Hz) the 
variables on se*(2) and by (111,112) the variables on R 2 , the product Poisson 
bracket is given by 



. df dg df dg 

{j,gh =mi « — p ~^—^r~ 

\Ofi2dfi3 d(x 3 dfX2 



V-2 



df_dg_ 
dfi 3 dfix 



df dg 
dfii d^ 3 



df_dg_ _ df dg 
du2 dux dux du2 
The Casimir function for this Poisson bracket is given by 

C(/xi,/z 2 ,M3) = 



2 



It is easy to verify that the linear transformation <E> : R 5 — * se*(2) x R 2 given 

by 

$(xi,X2,X 3 ,X4,,X 5 = (X1,X 2 ,X 3 ,X4,£X3 +X 5 ). 

is a Poisson diffeomorphism between (R 5 , {-, -}i) and (se*(2) x R 2 , {•, -} 2 ). 



In the new variables the system (2T) becomes 



/i 2 = nin 3 - e^iu 2 +£ 2 /iiM3 

/i 3 = -A*iM2 ( 2 - 2 ) 

Ui = -U 2 + £^3 
U 2 = Ul, 

which is an Hamiltonian system with respect to the Poisson bracket {•, -} 2 and 
the Hamiltonian function is given by 

H £ {m,fi 2 ,[i 3 ,ui,u 2 ) = -(ul + 2n 2 2 + i4 + ul + u 2 2 -2e^u 2 + e 2 yl). (2.3) 

3 Chaos by the Melnikov method. 



We will prove the occurrence of chaotic motion in the system (2.2) by showing 
the existence of transverse heteroclinic orbits. Melnikov Q gave an effective 
method to prove the existence of transverse heteroclinic (homoclinic) orbits in 
the Poincare map for a perturbed one-degree of freedom Hamiltonian system by 
measuring the "distance" between the stable and unstable manifolds associated 
with the saddle points. This method was generalized by Holmes and Marsden 
to the case of perturbed two-degree of freedom Hamiltonian systems when the 
phase space is a product of the dual of a Lie algebra and a set of action-angle 
variables. We briefly recall below this result; see Holmes and Marsden jl| for 
proofs. 

The setting is the following. The phase space is the product of the dual of 
a Lie algebra g and R 2 . The Hamiltonian has the form 

H E (n, 9, 1) = F(pi) + G(I) + eH\fi, 6, 1) + 0(e 2 ) 

= H a (fi,I)+eH 1 (n,0,I)+O{e 2 ) (3.1) 

where /i — (fi%, fi m ) £ $* and (0,1) are coordinates on R 2 , with 9 a 2ir- 
periodic variable. It is also assumed that the Lie-Poisson system whose Hamil- 
tonian is F has a heteroclinic (or homoclinic) orbit n(t) e g*. The oscillator 
frequency 

is assumed to be positive. The result is the following: 
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Theorem 3.1 (Holmes-Marsden) Suppose fi(t) is a heteroclinic (or homo- 
clinic) orbit for the Lie-Poisson system whose Hamiltonian is F, which lies in a 

two dimensional coadjoint orbit in g*. Let h — F(fi) be the energy of the hetero- 
clinic orbit and let h > h and 1° = G _1 (/i — h) be constants. Let {F, H 1 }^, 9°) 
denote the Lie-Poisson bracket of F(p) and H x (fi,n(l°)t + 6°,l°) evaluated at 
[i(t). Let 



oo 



and assume M(9°) has simple zeros. Then for e > sufficiently small, the 
Hamiltonian system l \3.f\ ) contains transverse heteroclinic orbits and hence Smale 
horseshoes on the energy surface H e = h. 

Now we will prove that the system (f^ ) verifies the conditions of the above 
theorem. The unperturbed system of ( |2.2| ) on se*(2) has unstable critical points 
(0, ±M, 0) lying on the 2-dimensional coadjoint orbit given by the cylinder 

{(Mi,M2, Ms) G se*(2) : M ? + £ = M 2 }. 

On this coadjoint orbit we have the heteroclinic orbits given by 

jUi(t) = ±Msech(Mt) 
fi 2 {t) = ±M tanh(Mi) 
fx 3 (t) = ±Msech(Mt) 

that link the unstable critical points (0, M , 0) and (0, —M, 0). 

On M 2 the unperturbed system of ( |2.2| ) is completely integrable and in action- 
angle coordinates (L,9), u\ — \f2Icos9, u 2 — \f2Isii\9, takes the form 

J = 
9= 1. 



Now writing the Hamiltonian (2.3) in the form 

H £ (^,fi 2 , fi 3 , 1,9) = - (nl + 2/x| + nl) + I - e^ 3 V2I sm0 + 0{e 2 ). 

wc arc in the setting of the Holmcs-Marsdcn theorem and we can write the 
heteroclinic orbits for an energy level 

H° = h = M 2 + k, 

where M 2 = /i 2 + [i 2 , and k is a constant, as 

/ui(t) = ±M sech(Mt) 
H2{t) = ±Mtanh(Mt) 
H 3 (t) = ±Msech(Mt) 
I = k 
9 = t + 9°. 
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The Melnikov function is: 



/CO 
V2kM 2 sech(Mi) tanh(Mi) sin(t + 6»°)dt 
-OO 

= -V2kM 2 (^j sinh(Aft) cosh" 2 (Aft) sin(t)d/j cos 6>° 



= — TrV^fcsech ( — ^— ] 
\2M) 



cos 6»°, 



which has simple zeros as a function of 8° and therefore the Hamiltonian system 
( fj!^ ) for e > sufficiently small, has transverse heteroclinic orbits and hence 
Smale horseshoes in a suitably chosen cross section of the constant energy surface 
with k > 0. 
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